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Abstract 

Here is studied the presence of absorption components shifted to the violet or the red side of the main 
spectral line (satellite, or discrete absorption components, i.e. SACs or DACs), in Mg II resonance lines' 
regions in Be stars and their kinematical characteristics. Namely our objective is to check if exists a 
common physical structure for the atmospheric regions creating SACs or DACs of the Mg II resonance 
lines. In order to do this, a statistical study of the Mg II AA 2792.523, 2802.698 A lines in the spectra 
of 64 Be stars of all spectral subtypes and luminosity classes is performed. We found that the absorption 
atmospherical regions where the Mg II resonance lines originated may be formed of several independent 
density layers of matter which rotate with different velocities. It is attempted also to separate SACs and 
DACs according to low or high radial velocity. The emission lines were detected only in the earliest and 
latest spectral subtypes. 
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1. Introduction 

The Mg II resonance lines have a peculiar profile in the 
Be stellar spectra which indicates multicomponent nature 
of their origin region. Many researchers have observed 
the existence of absorption components shifted to the vi- 
olet or the red side of the main spectral line (Underbill 
1970; Marlborough, Snow, & Slettebak 1978; Dachs 
1980; Doazan 1982; Danezis 1984; Danezis 1987; Danezis, 
Theodossiou, & Laskarides 1991; Sahade, Brandi, & 
Fontcla 1984; Sahade & Brandi 1985; Hutsemckcrs 1985; 
Doazan et al. 1991; Laskarides, Danezis, & Theodossiou 



1992; Cidale 1998; Lyratzi et al. 2003). These com- 
ponents, so called Discrete Absorption Components - 
DACs (Bates & Halliwell 1986) or Satellite Absorption 
Components - SACs (Danezis et al. 2003; Lyratzi & 
Danezis 2004), probably originate in separated regions 
which have different rotational and radial velocities. 
Especially in the case of very narrow DACs and SACs, 
they cannot be photosphcric, but rather they have cir- 
cumstellar or interstellar origin (Slettebak & Snow 1978). 
For example, Kondo, Morgan, & Modisette (1976) found 
that the shell absorption gets stronger from intermedi- 
ate to the late B stars and suggested that "it might be 
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due to the rising temperature in the gaseous sheU which 
converts Mg II to Mg III and to the weakening of the 
outward-driving mechanisms of the atmosphere" . In any 
case, the whole feature of the Mg II resonance lines is 
not the result of a uniform atmospherical region, but the 
components are created in different regions, which rotate 
and move radially with different velocities. As de Jager et 
al. (1979) proposed in their study of 33 stellar spectra of 
all the spectral types, in the late B supergiants variable 
mass loss occurs, due to "occasional stellar "puffs" super- 
posed on a more or less regular wind" . They proposed 
that "there are concentrations of low-ionization species in 
the stellar wind as a result of the occurrence of significant 
density variations" . Also, in order to explain the complex 
profiles of the Mg II resonance lines, Cidale (1998) pro- 
posed that the Be stellar atmospheres are composed from 
a classical photosphere, an extending high temperature 
chromosphere and a cool envelope. There is a question 
about contribution of different atmospherical layers, es- 
pecially in the case of SACs and/or DACs phenomena, in 
the construction of the Mg II resonance lines' profiles. 

The aim of this work is to statistically investigate the 
presence of SACs and/or DACs in Mg II resonance lines' 
regions in Be stars and their kinematical characteristics. 
Also, we would like to conclude about the limits of the 
rotational and radial velocities {Vrot, Vrad), as well as to 
check whether there exists a common physical structure 
for the atmospherical regions which create the Satellite 
Absorption Components (SACs) of the Mg II resonance 
lines in the spectra of all the Be stars. To do that, a statis- 
tical study of the UV Mg II resonance lines AA 2795.523, 
2802.698 A in the spectra of 64 Be stars (all the spec- 
tral subtypes and luminosity classes) has been performed. 
The study is based on the model proposed by Danczis et 
al. (2003) and Lyratzi & Danezis (2004). 

We decided to study the Mg II resonance lines, as they 
are characteristic of the cool envelope in Be stellar atmo- 
spheres, they are very intense features in the spectra of 
Be stars and they mostly present a complex and peculiar 
structure. Besides, since they are resonance lines, they 
give us a possibility to test the validity of the proposed 
model, as we have to adhere to all the necessary phys- 
ical criteria and techniques (Danezis et al. 2003). Our 
purpose here was to study spectral lines which are cre- 
ated in cool regions (Mg II) in Be stellar atmospheres. 
Danezis (1987); Sahade, Brandi, & Fontela (1984); Sahade 
& Brandi (1985) detected multistructure in the regions 
where Fe II lines (I.P.=7.870 cV) are created in the spec- 
tra of Be stars which are characterized as iron stars. One 
of our purposes was to investigate whether the multistruc- 
ture appears only in the Fe II spectral lines or in other 
spectral lines also with similar Ionization Potential, as Mg 
II (I.P.=7.646 eV). This is another reason why we chose 
to study the Mg II resonance lines and not some other 
resonance lines of the cool envelope, as N II (I.P.=14.490 
eV), C II (I.P.=11.260 eV), Si II (I.P.=8.110 cV) or Al II 
(I.P.:=5.986 cV). The study of Mg II resonance lines, gives 
us information not only about the cool envelope, but also 
for the multistructure of another ion that lies in the same 



region where the Fe II spectral lines are created. 

In §2 we describe the method of analysis, in §3 the ob- 
servational data, in §4 the results and their discussion is 
given and in §5 we give our conclusions. 

2. Method of spectral line analysis 

Before starting with the description of the method used 
in our study, let us explain the differences between DACs 
and SACs. The DACs arc components of a spectral line 
of a specific ion, shifted at different AA from the tran- 
sition wavelength of a line, because they are created in 
different density regions which rotate and move radially 
with different velocities (Lyratzi & Danezis 2004). The 
DACs are discrete lines, easily observed, in the spectra of 
some Be stars of luminosity class III, in the case of the 
Mg II doublet. However, if the layers that give rise to 
such lines rotate with quite large velocities and move ra- 
dially with small velocities, then the produced lines are 
quite broadened and slightly shifted. As a result, they are 
blended among themselves as well as with the main spec- 
tral component and thus they do not appear as discrete, 
consequently they cannot be resolved. In such a case the 
name Discrete Absorption Component is inappropriate. 
Besides, as Peton (1974) first pointed out, these compo- 
nents appear as "satellites" in the violet or in the red side 
of the main spectral line's component as a function of the 
time or the phase in the case of a binary system. For 
these two reasons and in order to include all these lines 
that are components of the line profile either they are 
discrete or not to a unique name, it has been proposed 
(Lyratzi & Danczis 2004) that they should be named 
Satellite Absorption Components (SACs) as a general ex- 
pression, because the Discrete Absorption Components 
arc "Satellites" of a main spectral line, while the Satellite 
Absorption Components are not always discrete and can- 
not be easily resolved. 

An additional peculiar phenomenon is that not all the 
lines of a specific ion have DACs. The DACs phenomenon 
is not a general one, but it is present in the case of some 
lines which have low excitation potential. For example, 
while the Mg II resonance lines at 2795.523, 2802.698 A 
(multiplet 1) have DACs in the spectra of some Be stars, 
their subordinate hues at 2790.768, 2797.989 A (multiplet 
3) do not present the same phenomenon. Therefore, we 
decided to study the Mg II resonance lines, as they are 
the only Mg II doublet having DACs. 

In order to try to obtain a qualitative picture we will use 
as an approximation the spherical geometry. According to 
it, there are two possibilities (see Fig. 1). 
A) The region which creates the SACs or DACs may be 
an envelope around and near the rapidly rotating star. In 
such a case, the spherical model could give only a rough 
approximation but allow however to obtain some quali- 
tative results. We note that such assumption is valid for 
some earlier models like e.g. Doazan and Thomas (1982) 
which can not explain the observed free-free emission or 
polarization, but we arc interested here only in kinematic 
characteristics. In this case the calculated values of the 
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Fig. 1. Density regions which create the observed SACs or 
DACs in the stellar spectra. 

radial velocity correspond to the component of the expan- 
sion or contraction velocity of the shell, which is projected 
to the observational axis. 

B) The region which creates the SACs or DACs may be 
an independent density region (blob) . which is spherically 
symmetric around its own center. Such density regions 
have been observed around stars that eject mass. In 
(Fig. 2) we sec this phenomenon in WR 104, observed 
by Tuthill, Monnier & Danchi (1999). Such a region may 
be the one, which creates the Mg II resonance lines and 
which lies in the cool envelope of the stellar atmosphere. 
Such a blob may have three different motions: a) it may 
rotate around the star, b) it may expand or contract and 
c) it may move radially. This means that the calculated 
values of the radial velocity consist of three different com- 
ponents: a) the component of the rotational velocity of 
the blob around the star, projected to the observational 
axis, b) the component of the expansion or contraction 
velocity of the blob, projected to the observational axis 
and c) the component of the velocity of the blob's radial 
motion, projected to the line of sight. 

In principle, as it is known, the star ejects mass with a 
specific radial velocity. The stream of matter is twisted, 
forming density regions such as interaction of fast and 
slow wind components, corrotating interaction regions 
(CIRs), structures due to magnetic fields or spiral streams 
as a result of the stellar rotation (Underbill & Fahey 
1984; MuUan 1984a; MuUan 1984b; MuUan 1986; Prinja 
& Howarth 1988; Cranmer & Owocki 1996; Fullerton et 
al. 1997; Kaper et al. 1996; Kaper et al. 1997; Kaper et al. 
1999; Cranmer, Smith, & Robinson 2000). Consequently, 
hydrodynamic and magnetic forces act as centripetal 
forces, resulting that the outward moving matter twists 
and moves around the star (see Fig. 2) . This motion of 
the matter is responsible for the formation of high den- 
sity regions (shells, blobs, puffs, spiral streams), which, 
either they are spherically symmetric with respect to the 
star or with respect to their own center (Danezis et al. 
2003; Lyratzi & Danezis 2004) (see Fig. 1). 

In order to study the physical structure and the exis- 
tence of SACs phenomena in the regions where these lines 
are created we used the model proposed by Danezis et al. 
(2003) and Lyratzi & Danezis (2004). This model allows 
us to calculate the rotational (K-ot) and radial velocities 
^rad) of independent density layers of matter in these 
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Fig. 2. Due to the rotation of the star the ejected matter 
forms spiral streams which produce blobs that create the ob- 
served SACs or DACs in the stellar spectra. This figure is 
taken from Tuthill, Monnier & Danchi (1999) 

regions, as well as the optical depth (r) and the column 
density (n). If the considered SACs or DACs originate 
from "puts" or "blobs" created by stellar winds in a cool 
extended envelope, the real shape of the envelope is not 
crucial. If the considered features are created in the en- 
velope's layers we also could obtain some qualitative con- 
clusions, since we are interested only in kinematical prop- 
erties here. Moreover, there is no a sophisticated or rough 
model of Be star which could fit the observed profiles con- 
sidering their complex structure so that there is no results 
of similar investigations for comparisons. 

Let us consider that the area of gas, which creates a spe- 
cific spectral line, consists of i independent absorbing den- 
sity regions followed by j independent regions that both 
absorb and emit and an outer absorbing region. Such a 
structure produces DACs or SACs in the observed spectra 
(Danezis et al. 2003) and the final line function which can 
describe the complex line profiles of the observed spectral 
lines is (Danezis et al. 2003; Lyratzi & Danezis 2004): 

where F\ is the observed flux, F\q is the initial flux, S\ej 
are the source functions of each emitting density region, 
^ is the optical depth in the center of the line and L^, 
Lej, i^g are the distribution functions of the absorption 
coefhcients /cai, A:Aej, ^Agi respectively. Each L depends 
on the values of the rotational and the radial velocity of 
the density region, which forms each component of the 
spectral line (K-ot, K-ad) (Danezis et al. 2003; Lyratzi & 
Danezis 2004). The product of L and ^ is the optical 
depth of each region. 

This function docs not depend on the geometry of the 
regions creating the observed feature. The considered ge- 
ometry of the regions is taken into account in order to de- 
fine the distribution function L. This means that L may 
represent any distribution which considers certain geom- 
etry (see Appendix A), without changing anything in F\. 

Each component of the spectral line, which is formed by 
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the density region of matter, must be accurately repro- 
duced by the function e^-^'^' by applying the appropriate 
values of Vroti, Vradi and ^i. Using the best model's fit 
for a complex spectral line, we can calculate the apparent 
radial and rotational velocity {Vradi , Vroti) and the optical 
depth (^i) in the center of the line of the region in which 
the main spectral line and its SACs are created. 

In the case that we want to consider that some other 
physical parameters are responsible for the line broaden- 
ing and not the rotation of the region which produces 
the studied spectral lines, we may replace the exponential 
e~^? with another classical distribution. 

In the case of emission lines, each emission compo- 
nent, which is formed by the j*'* density region of matter, 
must be accurately reproduced by the function Sxeji^ — 
6"^=^^=^), by applying the appropriate values of Vroti, 
Vradi, ^nd S. Usiug the best model's fit for a complex 
spectral line, we can calculate the apparent radial and ro- 
tational velocity {Vradi, Vroti), thc optical depth (^i) in the 
center of the line and the source function S of the region 
in which the emission component is created. 

The proposed model is relatively simple, aiming to de- 
scribe the regions where the spectral lines which present 
SACs are created. With this model we study the regions 
of a specific ion which creates a specific spectral line. 

This model presupposes that the main reason of the 
line broadening is the rotation of the region that gives 
rise to the spectral line (Doazan 1982). We can accept 
this assumption when we deal with the inner layers to the 
post-coronal regions. Thus, for these atmospheric layers 
the model gives satisfactory results. 

3. Observational data and fitting procedure 

The data we used are the Mg II resonance lines of 64 Be 
stars taken from the lUE Archive Search database^. The 
stellar spectra were observed with lUE satellite using the 
Long Wavelength range Prime and Redundant cameras 
(LWP, LWR) at high resolution (0.1 to 0.3 A). In Table 

I we give the list of stars, their spectral type and the type 
of camera used during observations. 

Our first step is to identify the spectral lines in the stud- 
ied wavelength range, in order to find out which lines may 
be blended with the Mg II doublet and, thus, may con- 
tribute to the observed features. The identification has 
been made by NIST Atomic Spectra Database^, as well 
as the catalogues of Moore (1968) and Kelly (1979). In 
this specific spectral range the adjacent features of the Mg 

II profiles are intense, but they are away from the Mg II 
spectral features, so that, in spite of their important influ- 
ence on the wings, their much smaller influence in central 
parts is not important for our discussion. Moreover, as we 
deal with resonance lines, we know that if one line of the 
doublet is well fitted, we should apply the same param- 
eters to the other one, even if the fit is not so good. In 
this case the unfitted regions correspond to blends. Also, 



the level of continuum is calculated for the whole spec- 
trum taken by lUE. This means that its suppression by 
the blends does not affect our calculations. 

The procedure that we followed for the decomposition 
of the lines is described in Danezis et al. (2003). The way, 
the criteria and discussion about the model we use to fit 
the observed spectral lines are presented in Danezis et al. 
(2003) and they are explained in more detail in Lyratzi 
& Danezis (2004) and an application of the method is 
given in Popovic ct al. (2004). This means that we tried 
to fit the observed profiles of the Mg II resonance lines 
with the less possible components. At first we tried to fit 
the observed spectral feature with one component. If this 
had not been possible, we added one more component and 
then another, until we accomplished the best fit with the 
less possible components. This means that we fitted the 
observed Mg II features with one to four components, for 
different stars. Our results are shown in Tables 2 and 3. 



http: / /archive. stsci.edu/cgi-bin/iuc 

http: / /physics. nist.gov/cgi-bin/AtData/lincs_form 
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Table 1. The list of Be stars with spectral type (columns 2, 5) and the type of the camera used during the observations (columns 
3, 6). 
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* Double system ** Triple system 

References: (1) Morgan, Code & Whitford 1955; (2) Lesh 1968; (3) Sahade, Brandi & Fontela 1984; (4) Herbig & 
Spalding 1955; (5) Guetter 1968; (6) Houk N. & Smith-Moore 1988; (7) Cowley 1972; (8) Osawa 1959; (9) Cowley et 
al. 1969; (10) Slettcbak 1954; (11) Houziaux & Andrillat 1976. 
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Table 2. The kincmatical parameters for absorption regions for considered stars: the rotational (Vrot) and the radial (Vrad) 
velocities in km s^^. 



Star 


Spectral Type 


Vrotl 


Vradl 


Vrot2 


Vrad2 
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Vroti 
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28 


-0.4 


















HD 193237 


B2 


pe 


30 


-18.6 


45 


-18.1 


86 


(DAC) 


-200.6 


115 


(DAC) - 


-200.1 


HD 45910 


B2 


HI e 


24 


13.0 


42 


11.0 


60 


(DAC) 


-229.5 


170 


(DAC) ^ 


-182.0 


HD 37202 


B2 


IV p 


29 


-18.0 


48 


-21.0 








125 




-33.0 


HD 36576 


B2 


IV - V e 


29 


-43.0 






80 




-21.0 








HD 212076 


B2 


IV - V e 


20 


-1.0 


34 


-1.5 














HD 32991 


B2 


V e 


19 


-15.4 


38 


-21.9 














HD 58050 


B2 


V e 


25 


-35.9 


49 


-28.4 














HD 164284 


B2 


V c 


28 


13.3 






63 




33.8 








HD 41335 


B2 


V ne 


22 


-9.5 


38 


-1.5 


80 




27.0 


100 




40.0 


HD 52721 


B2 


V ne 


18 


12.3 


34 


11.3 


65 




31.3 








HD 58343 


B2 


V ne 


21 


-1.8 


50 


-5.2 














HD 148184 


B2 


V ne 


17 


19.9 


30 


20.9 














HD 202904 


B2 


V e 


20 


4.0 


31 


4.5 














HD 65079 


B2 


V ne 


29 


8.0 






74 




-12.5 








HD 28497 


B2 


V ne 


8 


-45.0 


37 


-34.5 


60 




-35.0 








HD 45995 


B2 


V nne 


29 


-17.0 


45 


-20.0 


69 




-24.0 








HD 10516 


B2 


V pe 


22 


-1.0 


37 


-2.3 


50 




-5.8 








HD 187567 


B2. 


,5 IV e 


19 


23.6 


32 


24.6 


50 




26.1 








HD 32343 


B2. 


.5 V e 


20 


8.4 






78 




17.9 








HD 65875 


B2, 


,5 V e 


17 


14.0 


35 


18.0 














HD 191610 


B2, 


,5 V e 


20 


-45.6 


45 


-53.6 














HD 60855 


B2/B3 V 






45 


-14.6 


71 




22.9 


139 




-29.6 


HD 37490 


B3 


HI e 


20 


-5.3 


35 


-8.3 


54 




-18.8 


98 




-10.3 


HD 50820 


B3 


IV e+... 


24 


-13.0 


50 


-9.5 














HD 25940 


B3 


V e 


22 


0.7 


42 


-2.3 














HD 45725 


B3 


V e 


18 


-9.5 


34 


-10.5 


65 




-13.0 








HD 183362 


B3 


V e 


24 


24.2 


34 


24.2 


68 




19.7 








HD 208057 


B3 


V e 


23 


12.3 






75 




9.8 








HD 205637 


B3 


V:p 


22 


17.5 


39 


15.5 








115 




9.0 


HD 217543 


B3 


V pe 


26 


16.5 


51 


15.5 








130 




23.5 


HD 217050 


B4 


HI pe 


35 


14.3 


54 


17.8 


84 




32.3 








HD 89884 


B5 


HI 






33 


-9.5 


60 




1.5 








HD 22192 


B5 


V e 


16 


-2.0 


32 


-1.0 


52 




-7.0 








HD 23302 


B6 


HI e 


21 


-9.9 






86 




-8.9 








HD 45542 


B6 


HI e 


15 


-24.4 


33 


-24.9 


58 




-15.4 








HD 109387 


B6 


HI pe 


18 


-11.4 






66 




-28.4 








HD 23480 


B6 


IV e 


20 


-7.7 


32 


0.3 














HD 217891 


B6 


V e 


18 


-0.8 


















HD 138749 


B6 


V nne 


27 


22.5 


47 


14.5 














HD 23630 


B7 


HI 


20 


-9.0 


45 


-9.0 














HD 209409 


B7 


IV e 


21 


-9.0 


38 


-9.0 


83 




-5.0 








HD 6811 


B7 


V e 


19 


3.3 


35 


1.8 








169 




12.8 


HD 192044 


B7 


V e 


25 


19.5 


52 


8.5 














HD 210129 


B7 


V ne 


21 


62.6 


52 


43.1 














HD 142983 


B8 


la/ lab 


21 


-7.0 


39 


8.0 


60 




11.5 


95 




16.5 


HD 29866 


B8 


IV ne 


26 


-42.0 






85 




41.5 








HD 47054 


B8 


V e 


24 


-30.0 


42 


-31.5 
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Table 2. (Continued.) 



HD 


183914 


B8 V e 


17 


16.8 


37 


20.8 










HD 


50138 


B8 V e 


27 


-36.0 














HD 


58715 


B8 V var 


22 


-21.5 


44 


-9.0 










HD 


23552 


B8 V ne 






33 


24.7 










HD 


185037 


B8 V ne 


26 


12.0 


45 


18.0 










HD 


199218 


B8 V nne 


17 


4.7 


35 


4.2 






95 


7.7 


HD 


91120 


B8/B9 IV/V 


23 


-13.5 


39 


-15.0 










HD 


142926 


B9 pe 


18 


-12.9 


50 


2.1 


67 


7.6 


165 


3.6 


HD 


144 


B9 HI e 


28 


-30.6 


53 


(DAC) -188.6 


80 


(DAC) -187.1 


178 


(DAC) -129.6 



Table 3. The same as in Table 2, but for the emission component where it is present. 



Star 


Spectral Type 


^rote 


^rade 


HD 


203374 


BO IV pe 


20 


72 


HD 


45910 


B2 HI e 


34 


81 


HD 


32991 


B2 V e 


57 


13 


HD 


164284 


B2 V e 


120 


-1 


HD 


148184 


B2 V ne 


50 


33 


HD 


45995 


B2 V nne 


46 


102 


HD 


65875 


B2.5 V e 


30 


53 


HD 


50820 


B3 IV e+... 


59 


-87 


HD 


217891 


B6 V e 


50 


6 


HD 


192044 


B7 V c 


94 


7 


HD 


210129 


B7 V ne 


97 


-29 


HD 


47054 


B8 V e 


81 





HD 


50138 


B8 V e 


53 


120 


HD 


199218 


B8 V nne 


131 


15 
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Table 4. Values of the confidence with which the accepted fit is better than the fit with the one less component. 



Star 


4-3 comps 


6-2 


comps 


2-1 comps 


Star 


4-3 comps 


6-2 


comf 


HD 


5394 




0, 


.9132 




HD 


25940 




0, 


.9972 


HD 


53367 


1.0000 








HD 


45725 




1, 


.0000 


HD 


203374 


0.9483 








HD 


183362 


0.9807 






HD 


206773 




1, 


.0000 




HD 


208057 




1, 


.0000 


HD 


200310 








0.9993 


HD 


205637 


1.0000 






HD 


212571 




0, 


.9826 




HD 


217543 


1.0000 






HD 


44458 


0.9981 








HD 


217050 


1.0000 






HD 


200120 








1.0000 


HD 


89884 




1, 


.0000 


HD 


193237 


1.0000 








HD 


22192 




1, 


.0000 


HD 


45910 


0.1609 








HD 


23302 




1, 


.0000 


HD 


37202 


1.0000 








HD 


45542 




0, 


.9998 


HD 


36576 




0, 


.9931 




HD 


109387 


0.9996 






HD 


212076 




1, 


.0000 




HD 


23480 




0, 


.9998 


HD 


32991 




0, 


.6737 




HD 


217891 








HD 


58050 




0, 


.9978 




HD 


138749 




0, 


.9936 


HD 


164284 




0, 


.7826 




HD 


23630 




0, 


.9986 


HD 


41335 


1.0000 








HD 


209409 




0, 


.9954 


HD 


52721 




0, 


.9767 




HD 


6811 


1.0000 






HD 


58343 




0, 


.9928 




HD 


192044 




0, 


.8070 


HD 


148184 




0, 


.9396 




HD 


210129 




0, 


.8927 


HD 


202904 




0, 


.9475 




HD 


142983 


1.0000 






HD 


65079 


1.0000 








HD 


29866 




1, 


.0000 


HD 


28497 


1.0000 








HD 


47054 




0, 


.9086 


HD 


45995 


0.1169 








HD 


183914 




0, 


.9975 


HD 


10516 


0.9981 








HD 


50138 








HD 


187567 




0, 


.9747 




HD 


58715 




0, 


.9189 


HD 


191610 




0, 


.9990 




HD 


23552 








HD 


32343 




0, 


.7793 




HD 


185037 




0, 


.9500 


HD 


65875 




1, 


.0000 




HD 


199218 




0, 


.9693 


HD 


60855 


0.9800 








HD 


91120 




0, 


.9778 


HD 


37490 


1.0000 








HD 


142926 


1.0000 






HD 


50820 




0, 


.6221 




HD 


144 


1.0000 







2-1 comps 



0.6758 
1.0000 



No. 
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In order to be certain that we have accompUshed the 
best fit, we performed an F-test, between the fit that we 
accept as the best and a fit with one component less. We 
present the results of the F-tcst in Table 4, where we give 
the values of the confidence with which the accepted fit 
is better than the fit with the one component less. We 
did not perform an F-test between the best fit and a fit 
with one component more, as this last one fit presented 
extreme differences with the observed spectral line pro- 
files or because the values of the measured parameters go 
against the classical theory for resonance lines (Danezis et 
al. 2003; Lyratzi & Danezis 2004). 

4. Results and discussion 

DACs phenomenon is quite common in O- and early B- 
type stars. However, we found that DACs phenomenon 
is also observed in late B-type stars (as in the case 
of HD 144). Moreover, many researchers (Underbill 
1975; Morgan, Kondo, & Modisettc 1977; Marlborough, 
Snow, & Slettcbak 1978; dc Jager et al. 1979; Doazan 
1982; Sahade, Brandi, & Fontela 1984; Sahade & Brandi 
1985; Hutsemekers 1985), even if they do not use the 
name DACs, observed the same phenomenon. We point 
out also, that in Underbill (1970); Marlborough, Snow, & 
Slcttebak (1978); Dachs (1980); Doazan (1982); Sahade, 
Brandi, & Fontela (1984); Sahade & Brandi (1985); 
Hutsemekers (1985); Doazan et al. (1991); Cidale (1998) 
only the known DACs phenomenon was investigated. 
However, in the present paper the similar phenomenon of 
SACs is introduced and we investigate whether this phe- 
nomenon is able to explain the complex structure of Mg II 
resonance lines in the stellar spectra of all the Be spectral 
subtypes. Our result is that the SACs phenomenon is able 
to explain this complex structure. 

Using the model described above we find the best fit 
for the Mg II resonance lines of 64 Be stars given in Table 
1. In Fig. 3, we give the best fit of the Mg II resonance 
lines for two stars (HD 45910 and HD 41335). The first 
one (HD 45910) presents DACs, where the decomposition 
of the three different components is easy, while the second 
one (HD 41335) presents SACs, where, it is hard to de- 
compose the three different components, without the con- 
venient model. So, one can see that the model describes 
conveniently the Mg II complex profiles and the complex 
structure of the regions where these lines are created. In 
the studied lUE spectra, the interstellar lines are system- 
atically shifted to the red for -|~99 ± 16 km s~^. We used 
the Hipparcos catalogues ^ and we applied the corrections 
for the systemic velocity of individual stars (Smith 2001) 
and for the orbital motion of the stars which are members 
of binary systems. Our results, shown in Tables 2 and 3 
as well as in Figs. 4-8, are correspondingly corrected. 

We should note that in our sample 15 binary and mul- 
tiple systems are present (see Table 1). Considering that 
the origin of DACs and SACs, in principle might be dif- 
ferent than in single stars, at the beginning we separately 
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Fig. 3. The Mg II resonance line profiles fitted with the 
model for HD 45910 and HD 41335. The thick lines present 
observations, and the thin lines the best fit. The DACs (HD 
45910) and SACs (HD 41335) are present below. 



http: / /vizier. u-strasbg.fr/viz-bin/VizieR-3 



10 



Lyratzi et al. 



Roratloual \"el*cltie* 



1 IM 



a ccHnpcHiwnt 


b cQmfKHieiit 


c 'Coiupcxnenl 


d condone nt 








































1 ■ 1 1 1 1 r r 1 1 1 1 



Bfl- k: B4 BS b: BJ B7 HI b:: B9 b: B3 BS 



Spectral sublj-pe 



(5 ffluhtp]* STSS 







1 






iO.d 




[)« 




■»« 




-100.« 



M ao 



Fig. 4. Mean values for each spectral subtype of the rota- 
tional velocities of all the SACs as a function of the spectral 
subtype. 
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Fig. 5. Mean values for each spectral subtype of the radial 
velocities of all the SACs as a function of the spectral subtype. 
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Fig. 6. Radial velocities of all the SACs as a function of the 
respective rotational velocities. 
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Fig. 7. Mean values for each spectral subtype of the rota- 
tional and radial velocities of the emission component as a 
function of the spectral subtype. 



Fig. 8. Radial velocities of the emission component as a func- 
tion of the respective rotational velocities. 



analyzed this kind of stars (presented with open circles 
in Figs. 4 - 6). However, we found that there is no 
systematic difference between single and binary /multiple 
stars. Consequently, in further analysis we will not treat 
binary/multiple systems in a different way. 

From the fit we obtained the rotational (Vroti ) and ra- 
dial velocities {Vradi ) for each Mg II resonance line origi- 
nating within the region. 

In Tables 2 and 3 we present the kinematical param- 
eters for the absorption Mg II resonance line forming re- 
gions (Tabic 2), as well as for the emission ones (Table 3). 
As one can see, not all the studied stars have an emission 
component, but only the ones listed in Table 3. One could 
assume that these emission components are the emission 
part of P Cygni profiles formed by scattering. In that 
case, their wavelength and width does not represent the 
radial and rotational velocity. This could explain the ma- 
jority of positive radial velocities since in such a case even 
in outflowing wind, lines will be widened and red shifted. 
However, two cases with large negative radial velocities in- 
dicate that the real picture may be more complicated. We 
know that the disk models, in many cases, may produce 
theoretically only the shape of line profiles, but they are 
not able to fit them. So, in such cases, we may consider 
that the observed profile results from a different mecha- 
nism. In Table 2 we give the list of stars, their spectral 
subtype (columns 1 and 2, respectively) and the values of 
the rotational velocities (columns 3, 5, 7 and 9) and the 
radial velocities (columns 4, 6, 8 and 10) of the respective 
components. Also, in Table 3, we present the kinemat- 
ical parameters of the emission components. In column 
1 we present the list of stars, in column 2 their spectral 
subtype and in columns 3 and 4 the values of the rota- 
tional and the radial velocity, respectively. Let us point 
out, here, that the calculated values correspond to the 
regions which create the SACs or DACs. Especially, the 
obtained rotational velocities correspond to the rotation 
of the region around itself and not around the star. 

In Figs. 4 and 5 we present, separately, the rotational 
velocities and the radial velocities, respectively, of all the 
SACs as a function of the spectral subtype. In Fig. 5 
the radial velocity's values which correspond to the DACs 
are clearly seen. In Fig. 6 we present the radial velocities 
of all the SACs as a function of the respective rotational 
velocities. In Fig. 7 we present the rotational and radial 
velocities' mean values of the two resonance lines of the 
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emission component as a function of the spectral subtype. 
Each point refers to the mean value extracted for each 
spectral subtype. Finally, the radial velocities of the emis- 
sion component as a function of the respective rotational 
velocities are given in Fig. 8. 

The values we have calculated lie within a small range 
and we can obtain only the mean values of radial and 
rotational velocities and their standard deviation. The 
points in the diagrams correspond to the mean values of 
the velocities for each spectral subtype and the error-bars 
which appear in some of the diagrams to the standard de- 
viation. This standard deviation is not only a statistical 
error but it includes also the possible variation of the incli- 
nation axis and the error of the spectral classification, as 
the spectral classification which was based in the optical 
range could not be appropriate for the UV range. This 
means that the error-bars which appear in the diagrams 
include the statistical error as well as the dispersion of the 
values due to the different axis inclination. 

The reproduction of the Mg II resonance lines AA 
2795.523, 2802.698 A using the model where SACs are 
present, suggests that the atmospherical regions, where 
the Mg II doublet is created, may be described in a 
unique way for all the studied Be stars. This result con- 
firms the suggestion given by de Jager et al. (1979) that 
in Be stellar atmospheres exists a concentration of low- 
ionization species in the stellar wind and it is due to the 
occurrence of significant density variations. This result 
is also in agreement with the results of Morgan, Kondo, 
& Modisette (1977), who proposed that there are "signif- 
icant absorption features" on the left side of each reso- 
nance line, which are attributed to "additional absorption 
within the stellar extended atmosphere". These "signifi- 
cant absorption features" can be the SACs that appear in 
the spectra of the early type stars. Danezis et al. (2003) 
and Lyratzi & Danezis (2004) suggested that the peculiar 
phenomena observed in the spectra of Oe and Be stars, 
such as the SACs, are due to independent density re- 
gions in the stellar environment. Such regions may be 
structures that cover all or a significant part of the stel- 
lar disk (shells, blobs, puffs, bubbles) (Underbill 1975; de 
Jager et al. 1979; Henrichs 1984; Underbill & Fahey 
1984; Bates & Halliwell 1986; Grady et al. 1987; Lamers, 
Faraggiana, & Burger 1980; Waldron, Klein & Altner 
1992; Cranmer & Owocki 1996; Rivinius et al. 1997; Kaper 
et al. 1996; Kaper et al. 1997; Kaper et al. 1999; Markova 
2000), interaction of fast and slow wind components, 
CIRs, structures due to magnetic fields or spiral streams 
as a result of the stellar rotation (Underbill & Fahey 
1984; MuUan 1984a; Mullan 1984b; MuUan 1986; Prinja 
& Howarth 1988; Cranmer & Owocki 1996; FuUerton et 
al. 1997; Kaper et al. 1996; Kaper et al. 1997; Kaper et 
al. 1999; Cranmer, Smith, & Robinson 2000). This is the 
common theory which explains the DACs phenomenon in 
early type stars. We have found DACs phenomena in early 
B-type stars (e.g HD 45910), as well as in late B-type stars 
(e.g. HD 144). 

As one can sec from Table 2, the SACs phenomenon ap- 
pears to be a classical one for the Be stars. In Fig. 3 and 
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in Table 2 one can see that all the studied stars present 
discernible (DACs) or indiscernible (SACs) components of 
the Mg II resonance lines. The indiscernible components 
appear in the spectra of all the stars of luminosity classes 
IV and V and most of the stars of luminosity class III. 
It is very interesting that the SACs are observed as dis- 
crete lines (DACs) in the spectra of the three stars HD 
193237 (B2 pe), HD 45910 (B2 III e) and HD 144 (B9 HI 
e), because they present quite different radial shifts. This 
means that the regions which create these lines move radi- 
ally with relatively large velocities producing lines shifted 
enough to be easily observed in the spectra. On the other 
hand, in the case of all the rest studied stars, the SACs of 
the Mg II resonance lines present similar radial velocities, 
resulting to the SACs being blended among themselves. In 
this case, we can distinguish these lines by the systematic 
differentiations of the rotational velocities. 

The decomposition of the observed profiles for the Mg II 
regions in Be stellar atmospheres confirms the existence of 
independent density regions, since by using such a struc- 
ture, we were able to reproduce the resonance lines of Mg 
II in all the studied stars. This decomposition is physi- 
cally meaningful as it enables us to detect kinematically 
different regions with different rotational and radial veloc- 
ities, as well as the optical depth and the column density, 
for each of the Mg II independent density regions, which 
produce DACs or SACs. 

Our analysis shows that regions where the considered 
Mg II resonance lines originate ("blobs" and "puts" cre- 
ated by winds or cool extended envelopes) may consist 
of more independent density layers of matter with differ- 
ent kinematical properties (one to four in the analyzed 
cases). Wc identified them by the decomposition of ob- 
served Mg II lines in the number of components which fit 
best the observed profile. Namely, depending on particu- 
lar star, we obtain that SACs or DACs may be divided in 
several rotational velocity groups (average values for rota- 
tional velocity groups found to be present at 64 considered 
stars are 22±5 km s~^, 41±7 km s~^, 69±11 km and 
130±31 km s^^). The corresponding radial velocities are 
near zero (-3.3±20.3 km - for the first density region; 
-3.6±20.6 km s~^ - for the second; -1.0±21.8 km s~^ - for 
the third one and +4.0±22.7 km s~^ for the fourth one). 
In the spectra of the stars HD 193237 (B2 pe), HD 45910 
(B2 HI e) and HD 144 (B9 HI e) the SACs appear as 
discrete components (DACs) and the radial velocities of 
the third and the fourth density region are -205.7±21.7 
km s^^ and -170.6±36.6 km s~^, respectively. In the case 
of the star HD 144, the second component appears also as 
a discrete component. The radial velocity of the region, 
which creates this component is -188.6 km s^^ (Figs. 4, 
5 and 6). The observed velocity dispersion may be due 
to the different values of the rotational axis inclination of 
the regions where the SACs are created. The results pre- 
sented above confirm that the Mg II doublet is more or 
less stable for a given spectral type as Gurzadyan (1975) 
suggested. We did not find any variation of the velocities 
in the Mg II regions with the luminosity class, except in 
the case of the peculiar stellar spectra, where the SACs 
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appear as discrete lines (DACs), while Kondo, Morgan, 
& Modisette (1976) proposed that, apart from the differ- 
ence among spectral subtypes, there is probably difference 
among luminosity classes too. 

We assume that independent density regions corre- 
sponding to particular components of considered Mg II 
lines (one to four components corresponding to one to 
four regions in the case of 64 stars analyzed here), lie all 
in the cool stellar envelope. Depending on the tempera- 
ture, different ions with different Ionization Potential are 
created in different regions at different distances from the 
star. This means that the spectral lines observed in the 
spectra of Be stars derive from specific atmospherical re- 
gions, different among themselves. The ions which are 
created very close to the star lie in regions which present 
spherical symmetry around the star (case A in Fig. 1). 
On the other hand, the ions which are created at long 
distance from the star lie in regions which present spher- 
ical symmetry around their own center and not around 
the star (case B in Fig. 1). As the ionization potential 
of the Mg II ions is I.P.=7.646 cV, the Mg II ions can be 
created only at great distance from the center of the star, 
meaning in the disk, where spherical symmetry around the 
star cannot be accepted. This means that the Mg II ions 
lie at regions which present topical or apparent spheri- 
cal symmetry (case B in Fig. 1). As a result, the Mg 
II spectral lines and their SACs/DACs may derive only 
from such regions (case B in Fig. 1) and not from clas- 
sically spherical regions around the star (case A in Fig. 
1). This kind of density regions (blobs) have been pro- 
posed by many researchers (Underbill 1975; de Jager et 
al. 1979; Lamers, Faraggiana, & Burger 1980; Henrichs 
1984; Undcrhill & Fahey 1984; MuUan 1984a; MuUan 
1984b; MuUan 1986; Bates & Halliwell 1986; Grady et 
al. 1987; Prinja & Howarth 1988; Waldron, Klein & 
Altner 1992; Cranmer & Owocki 1996; FuUerton et al. 
1997; Rivinius et al. 1997; Kaper et al. 1996; Kaper et al. 
1997; Kaper et al. 1999; Markova 2000; Cranmer, Smith, 
& Robinson 2000) and are detected in many other cases, 
as in active stars (WR 104, see Fig. 2) and many quasars, 
as we observe DACs/SACs in their UV spectra (Danezis 
et al. 2006). This means that the density regions are a 
common phenomenon, observed in different levels. The 
fact that we found one to four components is accidental. 
In principle, there could be more or less. Although the 
number of components is different in different stars, we 
may conclude that the Mg II resonance lines forming re- 
gions present a complex structure. 

Our proposition that SACs phenomenon is responsible 
for the structure of Mg II lines means that we expect the- 
oretically that the Mg II components have similar radial 
velocities, within the range of statistical error (~ 10 km/s). 
The problem was how to distinguish them. The common 
idea (Doazan 1982) is that the radial velocity of the kine- 
matically independent regions is a function of the distance 
from the rapidly rotating Be star. According to that, our 
first thought was to distinguish these regions according 
to their rotational velocities, which is confirmed by our 
calculations. Namely, the regions which create the Mg II 



lines have similar radial velocities and different rotational 
velocities. 

In our sample we cannot detect high radial velocities, 
whatever method we use (the proposed model or any other 
classical method). The radial velocities of all the SACs, in 
all the studied stars are about km/s and only in the case 
of the stars HD 193237 (B2 pe), HD 45910 (B2 III e) and 
HD 144 (B9 III e), where we observe DACs, we calculated 
radial velocities with values between -130 km/s and -230 
km/s. We can detect the same phenomenon in the case of 
Si IV in a sample of 68 Be stars (values of radial velocities 
between -116 km/s and +25 km/s) (Lyratzi et al. 2006), 
as well as in the case of Ha in 120 Be stars (values of ra- 
dial velocities around km/s) (Lyratzi et al. 2005). These 
results indicate that in the atmospherical layers from the 
photosphere (very broad components of Ha) to the cool 
envelope (Mg II resonance lines) we cannot detect very 
high radial velocities. However, as we expected, the pres- 
ence of DACs in three stars of our sample indicate that 
from the regions near to the star towards the ones away 
from the star, the radial velocity increases, but it does not 
reach high values, as happens in the case of Oe stars. In 
the case of specific Be stars that present DACs with high 
radial velocities in some of their spectra, e.g. 59 Cyg and 7 
Cas (Doazan et al. 1989; Telting & Kaper 1994), we should 
study each one of them as exception of the classical rule. 
This proposition is based on the fact that the spectral 
classification was made in the optical range and may not 
apply in the UV spectral range (Walborn 1971; Walborn 
& Panek 1984; Walborn & Nichols-Bohlin 1987) (see also 
the SIMBAD database: http://simbad.u-strasbg.fr/sim- 
fid.pl). This means that the spectral classification in the 
optical and the UV range are not always in accordance. 
As a result, some early Be stars could be late Oe stars and 
we know that it is a common phenomenon to observe high 
radial velocities of SACs/DACs in Oe stars. However, in 
the case of Be stars, we do not observe the same phe- 
nomenon. This difference in the behavior of density re- 
gions in Oe and Be stars is very interesting and requires 
further investigation. Finally, we should pay attention in 
the way the radial velocities are calculated. The classical 
method considered that the whole observed feature corre- 
sponds to only one spectral line, meaning that the radial 
velocity was calculated by the displacement of the deeper 
point of the observed feature. Considering the SACs idea, 
the observed feature consists of a number of spectral line 
components. As a result, the deeper point of the observed 
feature is only the result of the synthesis of all the SACs. 
In this case we should calculate the radial velocity of every 
one of these components. 

An emission component is present in the spectra of BO, 
B2, B2.5, B3, B6, B7, B8 and B9 type stars (Fig. 7). 
This means that the emission does not appear in the 
spectra of the middle spectral subtypes of the Be stars 
(Kondo, Modisette, & Woff 1975). The radial velocity of 
the emission component decreases as the rotational veloc- 
ity increases. The emission component presents positive 
or negative radial velocities. If one takes into account that 
negative radial velocities exist, which is not in agreement 
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Fig. 9. Mean values of the optical depth ^ in the center of 
the line for all kinematically separated components of each 
resonance line, as a function of the spectral subtype. 

with the assumption that all emission components are the 
emission part of P Cygni profiles formed by scattering, 
one can assume the following. The obtained velocities cor- 
respond to the regions where the emission component is 
created as e.g. strings, blobs, puffs, bubbles. This means 
that the emission region may approach or move away from 
the observer and its different position and motion around 
the star is responsible of whether this value is positive or 
negative. In Fig. 7 and 8 one can see that, as the rota- 
tional velocity increases, the radial velocity decreases, in 
contrast with the relation of the two velocities of the ab- 
sorption components. A problem with the emission com- 
ponent is that it is blended with absorption lines of other 
ions and thus it is difficult to evaluate the rotational and 
radial velocities. As a result the calculated values present 
greater statistical error than in the case of the absorption 
components. 

According to criteria described in Danezis ct al. (2003), 
the same components of the two resonance lines should 
have the same values of rotational and radial velocities 
and the ratio of the optical depth (^) should be the same 
as the respective ratio of relative intensities. In Fig. 9 
we present the mean values of the optical depth (^) in the 
center of the line for all kinematically separated compo- 
nents of each resonance line, as a function of the spectral 
subtype. As one can see, the value of ^ of the first kine- 
matical region is obviously higher than in the rest three. 
It tends to decrease from the first to the fourth kinemat- 
ical region. That, also, may indicate physically separated 
regions. 



region. We obtained the rotational and radial velocities 
which allow us to extract some general physical proper- 
ties for the Mg II regions of Be stars. Some interesting 
results inferred from the investigations are the following: 
(i) the proposed rotation model gives satisfactory results 
for the Mg II AA 2795.523, 2802.698 A resonance lines re- 
gion; (ii) the absorption atmospherical region where the 
Mg II resonance lines are created presents a complex struc- 
ture. It tends to be composed by more than one kinemat- 
ically independent regions (Only four stars present sim- 
ple structure). We found that the kinematically indepen- 
dent regions rotate with different velocities: 22 km s~^, 
41 km s~^, 69 km s^^, 130 km s^^. The respective radial 
velocities are near zero for all these regions. These calcu- 
lated values lead us to accept that the Mg II resonance 
lines of the Be stellar spectra present Satellite Absorption 
Components, (iii) the rotational velocities of the found in- 
dependent regions present a uniform fluctuation with the 
spectral subtype, (iv) the emission lines were detected 
in the earliest and latest spectral subtypes with the ma- 
jority of positive radial velocities with several exceptions, 
i.e. with radial velocities within the range to -87 km s~^ , 
and taking into account that the positive values, if lines 
are formed by scattering, are not radial velocities. 
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5. Conclusions 



We applied the method developed in Danezis et al. 
(2003) and Lyratzi & Danezis (2004) on the Mg II reso- 
nance lines of 64 Be stars in order to investigate the kine- 
matical properties of the Mg II resonance lines forming 
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Appendix 1. Calculation of the distribution functions L 

As we know the distribution function (i) of the absorption coefficient {k\) has the same form as the distribution 
function of each component of the spectral line. This means that we can replace the distribution function of the 
absorption coefficient Lj, with another expression of the distribution function of each component. 

It is also known that Be and Oe stars are rapid rotators. This means that we accept that the main reason of the line 
broadening is the rotation of the regions that produce each satellite component of the whole observed spectral feature. 
These rapidly rotating density regions may also present radial motion. For these two reasons we search another 
expression for the distribution function of the spectral line's components that has as parameters the rotational and 
radial velocities of the spherical density regions. 

For a spherical density region, we assume the following hypotheses: i) The natural broadening of the spectral lines 
follows the Lorentz distribution; ii)Lamberts' sinus low stands for each point of the spherical region; iii) The angular 
velocity of rotation is constant. 

In order to calculate the total radiation, we divide the spherical layer in very thin cylindrical surfaces which are 
perpendicular to the rotational axis. Lambert's low allows us to consider that the luminosity from each point on the 
sphere is the same. 

On the above cylindrical surfaces we also consider the surface dS. According to Lambert's low. when this surface 
rotates with an angular velocity w, its radiation intensity is: 



Ci is a constant and 7 is the Lorentzian full width at half maximum, which, in the case of the natural broadening, 
has the value 7 = IQ^Hz. 

When the surface dS does not rotate, the center of the formed spectral line has the observed wavelength Ao, which 
corresponds to a frequency i^q. Thus, we have: 
c 

ujQ = 2t:vq = 27r— . 

When the surface dS rotates with a rotational velocity Vrot, the center of the formed line has the wavelength 
and in this case we have uik = i^oi^ ~ zsiiup), where z = 

We also have cos 6* = cos a cos The angles a and ip are shown in Fig. 10. 



dl{uj) = Q{oj)dScoa0, 
where 8 is the angle between the vertical on dS and the line of sight and 



(Al) 
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Fig. 10. Elementary ring of the spherical density region. 



The surface dS can be written as dS = rdhdip, where r is the radius of the cylinder, dip is the angle under which 
the observer sees dS and dh is the height of dS. 
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Making the above substitutions in Eq. (Al), we have: 

Cirdhj cos a cos (pdip 

dl(uj) = 2 

[u - - zsi\iip)f + (^) 

Thus, the radiation intensity from the semi cyUnder is 

7r/2 

[ Cirdh"/ cosacosipdip 

= / 2 

~i/2 [^-^o(l-zsinv3)]2 + (2) 



■it/2 

ACirdh f cos ac? sin 



7 

-7r/2 



^,-,^)(l-zsinc,) 



2 

+ 1 



If we take that u) ~ ujq — x — sinos, we have 
72 72 

m 



1 

ACircosadh f dx 



- Cboil - zx)]^ + 1 



UJ 



Taking that y = uj — ujq{1 — zx), the above integral becomes 
AC ir COS adh f dy 

I{u)) 



). (A2) 



Finally, we have 

' 4Ci r cos ad/i \ / arctanp — a)o(l ^ ^)] ^ arctanp — a)o(l + z)] ' 
7 / V 

The above function describes the radiation intensity from a visual semi cylinder with radius r and height dh. 
Under the angle da is seen dh from the center of the spherical region. This cylinder rotates with a rotational velocity 
_ ^ constant angular velocity Cj. 

We consider the function: 

arctan[a) — a>o(l ^ 2)] — arctanp — a)o(l + -z)] 
F[uj,z) = z 

UqZ 

Wc study the limit of this function in the case that the density layer does not rotate, i.e. when ^ — ^ 0. In such a 
ise: 

arctanp — a)o(l — -z)! — arctanfa) — a)o(l + z)! 
lim_P(ti;,z) ~ Imi ; . 

We apply the De L' Hospital's low and we have: 

-^{arctanp-cjoll -z)]}- ■£{arctan[(:)-tJo(l + ^)]} 
hvaPiuj.z) = hva— ^ — 

1 1 



limP((I',z) — lim 



1 + P- (1.0(1 -^F l + p-il-oll + z) 



(j)-c;)o)2 + i' 



It is obvious that in the non-rotating case this form corresponds to the Lorentz's distribution for the naturally 
broadened spectral lines. 

In the case when the rotation broadening |Ai — A2I (or \uji — ^2!) is much larger than the natural broadening (the 
natural broadening of a spectral line is of an order of 10^'^ - lO^^A), the above function P{Cj,z) presents the form of 
one quadratic pulsation (see Fig. 11) 



This means that AXtotai = Ai — A2 = tzi§ and so 



For each iv with uji <uj <lj2 the relative shift is z = - — —. For the point oji we have "° = 2 so that uii = wo(l — 2). 
Likewise, for the point L02 we have 102 = t^o(l + z). But uj = 2itv = ^y^. Thus, Ai = and A2 = jr- with Ai > A2. 



l-z^ 

Amin ^ ^2 = Aq = Aq — Aq and Aniax = '^l = Aq H = Aq + Ag — 
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Fig. 11. Quadratic pulsation. 



We set p= Xojz^ and normalize to 1. In this way the function P{Cj,z) could be approximated with the function 
/(A) where: 

/(A) = ( I \^-M<P 
^ ' otherwise 



Now, we assume that the spherical density region rotates with equatorial velocity zq = — (Fig. 12) 



Obsen"er 



Fig. 12. Rotating spherical density region that produces SACs or DACs. 

The points of the circle with radius r rotate with a velocity Vrot and for r = i?cosa we take 
Vrot ^ Lor = ujRcosa. 

We set Vb = luR. Also uj ~ const. So: Vrot = ^Rcosa = Vbcosa and z = zocosa. We have also dh = Rda and Eq. 
(Al) becomes 



dI{LD) 



ACiR arctanp — i1}q{1 — zoCOsa)] — arctanp — a)o(l + ^ocoso;)] 



7 



ujqZo cos a 



The integral of this equation is 

7r/2 

x/~N 4Cii? f arctanfo) — a>o(l — ^ocosa) — arctanfa) — u;o(l + 2ocosa)l 
I[ujj = / ; ^cosaaa. 



7 



ujqZo cos a 



(A3) 



-7r/2 



When we take into account the function P(a),z) the above function becomes 

tt/2 



ACiR 

7 



j P [Co, zq cos a) cos ada. 



-7r/2 



We approximate P(a3, zpcosa) with /(A) for and we take the integral for the observation angle 6 <E [— 6*0, 6*0] from 
the equatorial plane. 
So we have that 



= Ii = / 1 • cos 
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h = 



ACxR 

1 



cos 9d9 : 



4Cii? 

7 



7 



- sm Wq . 



If we normaUze to 1 , we have 
h 



cos^ 6*0- 



sin 6*0 = yl^ 
For the angle we have 

AqZoCOS^o 



|A- Aol < po ^ 



1 



For a wavelength A or for a shift AA = |A — Ao| from the center of the spectral line, the absorbing (or emitting) 
regions are those with angular distance 9 from the equatorial plane, with \9\ < 9q. 
For the equatorial plane we have 

XqZq COS(?o 



AA: 



1 



cos^ 



From this equation we can calculate the angle 6*0 as follows: 



COSt/Q 



-Ao±v/A2+4AA2 



2AAzo 

As 9o is between — 7r/2 and 7r/2 we have cos9o > and so 



cosf^o = 



-Ao + VA^+4AA2 



2AAzo 

Thus, the distribution function Ii takes its final form: 



h = V 1 — cos-^ 9o if cos^^o = ""^zaa^T""" ^ ^ ^^"^ 
/i = 0, otherwise. 

It is obvious that the distribution fmiction Ii is a function of the wavelength (A). This means that Ii = /i(A). This 
distribution function has the same form as the distribution function of the absorption coefficient L and may replace 
it (in e~^^), when the main reason of the line broadening is the rotation. We name it Rotation distribution function. 

The spectral line profile, which is formed by a spherical density region, is reproduced by the function e^^^ by 
applying the appropriate value of the rotational velocity Vrot (from zq), the radial velocity Vrad (from ^^^^ = ^^^XTt^-' 
and the optical depth ^ in the center of the line. 
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